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Abstract 

We perform a rigorous analysis of the quasineutral limit for a hy- 
drodynamical model of a viscous plasma represented by the Navier 
Stokes Poisson system in 3 — D. We show that as A — the velocity 
field strongly converges towards an incompressible velocity vector 
field u and the density fluctuation — 1 weakly converges to zero. In 
general the limit velocity field cannot be expected to satisfy the incom- 
pressible Navier Stokes equation, indeed the presence of high frequency 
oscillations strongly affects the quadratic nonlinearities and we have to 
take care of self interacting wave packets. We shall provide a detailed 
mathematical description of the convergence process by using microlo- 
cal defect measures and by developing an explicit correctors analysis. 
Moreover we will be able to identify an explicit pseudo parabolic pde 
satisfied by the leading correctors terms. Our results include all the 
previous results in literature, in particular we show that the formal 
limit holds rigorously in the case of well prepared data. 
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1 Introduction and plan of the paper 



1.1 Introduction 

In this paper we perform a rigorous analysis of the so called quasineutral 
limit for a hydro dynamical model of a viscous plasma represented by the 
Navier Stokes Poisson system in 3 — D, namely 

dtp^ + div(p^n^) = 0, (1) 

dt{p^u^) + dw{p^u^(g,u^) + V{p^y = pAu^ + {u + p)Vdivu^ + p^VV^, (2) 

X^AV^ = p^ -I. (3) 

Let us denote by x G M^, t > 0, the space and time variable, p{x, t) the neg- 
ative charge density, m{x,t) = p{x,t)u{x,t) the current density, u{x,t) the 
velocity field, V{x,t) the electrostatic potential, p,,v the shear viscosity and 
hulk viscosity respectively. The parameter A is the so called Debye length 
(up to a constant factor). 

We show that as A — )• the velocity field strongly converges towards an in- 
compressible velocity vector field u and the density fiuctuation p"^ — 1 weakly 
converges to zero. In general the limit velocity field cannot be expected to 
satisfy the incompressible Navier Stokes equation, indeed the presence of 
high frequency oscillations strongly affects the quadratic nonlinearities and 
we have to take care of self interacting wave packets. In the paper we shall 
provide a detailed analysis of the convergence process by using microlocal 
defect measures and by developing an explicit correctors analysis. Moreover 
we will be able to identify an explicit pseudo parabolic equation satisfied by 
by the leading correctors terms. 

The previous system can be seen as the coupling of the compressible 
Navier Stokes equations ([T|) , ([2]) with a Poisson equation ^ , where in dimen- 
sionless units the coupling constant can be expressed in terms of a parameter 
A which represents the scaled Debye length, which is a characteristic physi- 
cal parameter related to the phenomenon of the so called "Debye shielding" , 
|14j . studied by Peter Debye in 1912. Any charged particle inside a plasma 
attracts other particles with opposite charge and repels those with the same 
charge, thereby creating a net cloud of opposite charges around itself, this 
cloud shields the particle's own charge from external view and then causes 
the particle's Coulomb field to fall off exponentially at large radii, rather 
than falling off as 1/r^. So the physical meaning of the Debye length A is 
the distance over which the usual Coulomb field is killed off exponentially 
by the polarization of the plasma. In terms of physical variables the Debye 
length can be expressed as 

A = Xn/L Xn = ./^^, (4) 
V e^no 
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where L is the macroscopic length scale, £q is the vacuum permittivity, ks 
the Boltzmann constant, T the average plasma temperature, e the absolute 
electron charge and uq the average plasma density. In many cases the Debye 
length is very small compared to the macroscopic length A/j << L and so it 
makes sense to consider the quasineutral limit A — t- of the system ([I])- ([3]). 
In this situation the particle density is constrained to be close to the back- 
ground density (equal to one in our case) of the oppositely charged particle. 
The limit A — )■ is called the quasineutral limit since the charge density 
almost vanishes identically. The velocity of the fluid then evolves accord- 
ing to the incompressible Navier Stokes flow. This type of limit has been 
studied by many authors. In the case of Euler Poisson system by Cordier 
and Grenier [^, Grenier [16], Cordier, Degond, Markowich and Schmeiser 
[3] , Loeper [2l] , Peng, Wang and Yong [25] , in the case of the Navier Stokes 
Poisson system by Wang [32] and Jiang and Wang [18] and in the contest 
of a combined quasineutral and relaxation time limit by Gasser and Mar- 
cati in [9l [TOl [TT] . This paper is still a mathematical theoretical approach 
to this complicate physical problem which however removes many regularity 
and smallness assumptions of various papers in the literature see for instance 
Wang [32] and Jiang and Wang [18] . In fact Wang [32] studied the quasineu- 
tral limit for the smooth solution with well-prepared initial data. Wang and 
Jiang [18] studied the combined quasineutral and inviscid limit of the com- 
pressible Navier- Stokes-Poisson system for weak solution and obtained the 
convergence of Navier- Stokes-Poisson system to the incompressible Euler 
equations with general initial data. Moreover in [18] the vanishing of vis- 
cosity coefficient was required in order to take the quasineutral limit and 
no convergence rate was derived therein. The authors in [6] investigated 
the quasineutral limit of the isentropic Navier-Stokes-Poisson system in the 
whole space and obtained the convergence of weak solution of the Navier- 
Stokes-Poisson system to the weak solution of the incompressible Navier- 
Stokes equations by means of dispersive estimates of Strichartz's type under 
the assumption that the Mach number is related to the Debye length. Ju, 
Li and Wang [19] studied the quasineutral limit of the isentropic Navier- 
Stokes-Poisson system both in the whole space and in the torus without the 
restriction on viscous coefficient with well prepared initial data. However 
there is no analysis for the quasineutral limit for the Navier Stokes Poisson 
system in the context of weak solutions and in the framework of general ill 
prepared initial data. The common feature of this kind of limits in the ill 
prepared data framework is the high plasma oscillations, namely the pres- 
ence of high frequency time oscillations along the acoustic waves. In these 
phenomena there are different behaviors of the various vector fields acting in 
our system. Particularly relevant us to understand the relationship between 
high frequency interacting waves, dispersive behavior and the different role 
of time and space oscillations. In our analysis the velocity fields both dis- 
perse and oscillates however the dispersion behavior dominates on the high 
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frequency time oscillations and Strichartz estimates are sufficient to pass 
into the limit of the convective term. The presence of quadratic terms on 
the electric field (e.g. p^\/V^) cannot be analyzed in the same way since 
the dispersive behavior no longer dominates on time hight frequency wave 
packets. In the general case these quadratic terms will not vanish in the 
limit as A — t- 0, unless we have well prepared initial data. 

1.2 Plan of the paper 

The structure of this paper, as well as the main ingredients of our approach 
to this limiting process can be summarized as follows. 

• In Section 2 we collect many needed mathematical tools, including 
notations, Strichartz estimates and microlocal defect measures. Then 
in Section 3 we set up our problem. 

• The following section 5 is devoted to obtain a priori estimates indepen- 
dent of A, namely standard energy bounds and dispersive estimates on 
the density fluctuation. The main idea here is based on the observation 
that the density fluctuation — 1 satisfies a Klein-Gordon equation, 
so acoustic waves analysis for the Navier Stokes Poisson system ([TJ- 
^ follows by reading the system as a dispersive equation and we will 
get uniform estimates in A by the use of he L^-type estimates due to 
Strichartz \13\ \20\ [28] . The particular type of Strichartz estimates for 
the Klein Gordon equation that we use here can be recovered from the 
seminal paper by Strichartz [25] (where he studied the homogenous 
equation) and by using Duhamel's principle. 

• In the previous sections we get sufficient bounds in order to study the 
limiting behaviour of the velocity vector field. Therefore in Section 
5 we analyze separately the limiting behaviour of the divergence free 
part and the gradient part of u^. Accordingly we obtain the strong 
convergence of the velocity field 

• The next stumbling block is to get enough compactness for the electric 
field in order to pass into the limit in the quadratic term XVV^ 
XV V. Since XVV^ is bounded in L^Ll we can define microlocal 
defect measure introduced by P. Gerard in |12] and by L.Tartar (H- 
measure) in [29] with correctors E'^ and E~ to handle time oscillations 
at frequency 1/A. An analogous use of the P. Gerard and L. Tartar 
ideas can be found in Y. Brenier and E.Grenier [2] and E. Grenier |15j . 
regarding the Vlasov Poisson system. This will be done in Section 6. 

• In Section 7 we will be able to prove our Main Theorem [TJ 

• As a final step, in Section 8, we show that in the case of smooth 
solutions for the system ([I])-(l3|) the class of correctors and E~ is 
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not empty and they satisfy a "pseudo parabolic" type equation, see 
the Main Theorem [2l 

2 Preliminaries 

For convenience of the reader we estabhsh some notations and recall some 
basic facts that will be useful in the sequel. 

2.1 Notations 

If F, G are functions we denote by -F < G the fact that there exists c G M 
such that F < G. Then, we will denote by 

a) P (M'^ X M+ ) the space of test function (M'' x ]R+ ) , by (M'' x R+ ) the 
space of Schwartz distributions and (•, •) the duality bracket between 
V and V 

b) W'''P{R'^) = {I- A)~tLP(M'^) and H''{R'^) = W'^'^iR'^) the nonhomo- 
geneous Sobolev spaces, for any 1 < p < oo and /c G M. VF'^'P(R'^) = 

{-Ay^LPiR"^) and H''{R'^) = W'''^{R'^) denote the homogeneous 
Sobolev spaces. The notations L^Ll and L^Wx'"^ will abbreviate re- 
spectively the spaces LP([0, T]; L9(M'^)), and LP{[0,T];W''^''{R'^)). 

c) L^iR'^) the Orhcz space defined as follows 

LP{R') = {/ G LUR'') I |/|X|;|<1 G L\R'^), |/|x|;|>i G L^(M'^)}, 

(5) 

see [T], [52] for more details. 

d) the space of bounded operators, /C(M^) the space of compact 
operators, 

e) if X, Y are Banach spaces, C{X, Y) is the space of bounded operators 

f) Q and P respectively the Leray's projectors Q on the space of gradients 
vector fields and P on the space of divergence - free vector fields. 
Namely 

Q = VA-^ div, P = I-Q. (6) 

It is well known that Q and P can be expressed in terms of Riesz 
multipliers, therefore they are bounded linear operators on every W^'P 
{1 < p < oo) space (see [27]). 
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Next we recall the basic notations concerning pseudo-differential operators 
and symbols to be used later on. We refer to [31] for details. Assuming 
p, 5 £ [0,1], m £ M, we denote the set of symbols satisfying 

for all a, /3, where (^) = (1 + ICP)"*^^^- Ii^ such case we say that the associated 
operator denoted by OP{p{x,S)) is given by 

P{x,D)f{x) = j p(x,e)-F/(Oe"«de := OP{p{x,i)) 

(whereJ-'/(^) = (27r)~" J f{x)e~^^^dx denotes the Fourier transform of the 
function /) belongs to OPS^^. If there are smooth symbols pm-j{x,£,), ho- 
mogeneous in ^ of degree m—j for |^| > 1, i.e. pm-j{x, r^) = r"^~^pm-j{x, ^) 
for r > 0, 1^1 > 1, and if 

P{x,0 ~ ^Pm~j{x,^) 

in the sense that 

N 

for all N, then we say p{x, ^) € S*"* and P{x, D) is polyhomogenous of order 
m. If is an open set in M^, we denote by ilj^^p{^, C{H)), respectively, 
'4'comp{^,f(^{H)) the space of polyhomogenous pseudo-differential operators 
of order m on with values in C{H), respectively IC{H) whose kernel is com- 
pactly supported in 17 x Jl, moreover we recall that if P G VSmpC^j '^(-^))' 
then its symbol p[x,S) is a linear application from ip^^pi^, C{H)) to 
C§°{S*n,C{H)), where S*Q = S"^-^ x n. 

Following P. Gerard we say that /x is the microlocal defect measure (or 
following L. Tartar the H-measure) for a bounded sequence Wk in if for 
any A G t/j'^^^^ploo , ]C{H)) one has (up to subsequences) 

lim {A{wk — w), {wk — w)) = / tr{a{x,£)ii{dxd£)). 

Js'-Q. 

where A = OP{a{x,£)). 
2.2 Technical tools 

2.2.1 Strichartz estimates for Klein Gordon equations 

Let us recall that if u; is a solution of the following Klein Gordon equation 
in the space [0, T] x 

-— A - j w{t, x) = F{t, x) 
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with Cauchy data 

w{0,-) = f, dtw{0,-)=g, 

where m > is the mass and < T < oo, then w satisfies the following 
Strichartz estimates, (see |28| ) 

II^^IIl?,, + \\dtw\\^,^-i., < 11/11^1/2 + + IIFIL?,,, 

where {q,p), are admissible pairs, namely they satisfy 

2(n + l) 2(n + 2) 2(n + 2) 2(n + l) 

^— ^<P<^— 7^ ^ -<P<- T^- 

n + 6 n + 4 n n — 1 

In particular in the case of d = 3, {q,p) are admissible if they satisfy 

4 10 10 4 

- <p< — —<„<-_ 

Moreover by choosing p = 4/3 and q = 4 and by a standard application of 
Duhamel's principle we have the following estimate 

\\w\\li^ + \\dtw\\^^^-iA < 11/11^1/2 + lbll^i/2 + \\F\\lIlI- (7) 

It is straightforward to observe that for any s > also this estimate hold 

MliW-^'* + \\9tw\\^4w-l-s,4 < 11/11^1/2-. + llff 11^-1/2-. + WHlIH-^- (8) 

(it is sufficient to apply the operator (/ — A)~^/^ to ([7])). 

2.2.2 Properties for pseudo-differential operators 

We recall here two fundamental tools necessary to work with pseudodiffer- 
ential operators (fore more details see [31], [3D], [12], ) 

Proposition 2.1. If A € OPS'^, then 

A:LUn,H)^L\n,H) 

is bounded. 

Proposition 2.2 (Generalized Rellich Theorem). If A G tlj^^^p{n,]C{H)) 
for some m < 0, then 

A:LUn,H)^L\n,H) 
is compact, i.e. if Wk ^ w weakly in L^, then \\Awk — w|| — )• strongly. 
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2.2.3 Convolution estimate and compactness theorems 

Here we state the following elementary lemma that will be used later on. 

Lemma 2.3. Let us consider a smoothing kernel j G C^(M'^), such that 
j > 0, J^djdx = 1, and let us define 

ja{x) = a-'^j (^^) . 
Then for any f £ H^{W^), one has 

\\f-f*ja\\LPiR'^)<Cpa~i-^~^)\\Vfh2^^,), (9) 

where 

pe[2,oo) ifd = 2, pe[2,6] if d = 3. 
Moreover the following Young type inequality hold 

\\f*ja\\LPiR'i) < Ca'"'^(^~^)||/||^-..,(Rd), (10) 

for any p,q £ [l,oo], q < p, s > 0, a £ (0, 1). 

We recall also the following compactness tool (see |26|). 

Theorem 2.4. Let he F C L'p{[Q,T];B), 1 < p < oo, B a Banach space. T 
is relatively compact in LP([0,T];S) for \ <P < oo, or in C([0,T];i?) for 
p = oo if and only if 

(i) <! J f{t)dt, f G is relatively compact in B , < ti < t2 < T, 
(ii) lim \\f{t + h) - f{t)\\LP([o,T-h];B) = uniformly for any f £T. 

3 Statement of the problem and Main Results 

3.1 Basic facts on the Navier Stokes Poisson System 

In order set up our problem we recall here some results concerning the 
existence theory for the Navier Stokes Poisson system ([I|)-(l3]). For simplicity 
we rewrite here again the compressible Navier Stokes equation coupled with 
the Poisson equation 

+ div(p^ti^) = 

< dtip^u^) + div(p^u^ u^) + V(p^)^ = fiAu^ + {i^ + fi)Vdivu^ + p^VV^ 
X'^AV^ = / - 1. 

(11) 
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To simplify our notation from now on we will set 

= ^^T^ = . = 

The system (jlip is endowed with the following initial conditions, 

pio = Po > 0, V\=o = V^, (ID) 

p^u\=Q = niQ, mo = on {x G I Pq{x) = 0}, 
,A|2 



/ (n%=o+^-^ + X'\Vo'\)dx<Co. 



The existence of global weak solutions for fixed A > for the system (jlip . 
has been proved in the case of a bounded domain in and in the case of 
the whole domain in [7] and [8]. We summarize this existence result in the 
following theorem. 

Theorem 3.1. Assume (ID), and let 7 > 3/2, then there exists a global 
weak solution {p^,u^,V^) to (dH) such that - 1 e L°°((0, T); L](]R3))^ 
e L°°((0,r);L2(]R3))^ g L2((0,r); W^1'2(]R3)). Furthermore 

• The energy inequality holds for almost every t > 0, 

Z^ + vr^A^lVy^l^) dx 

+ [ f (p\Vu^\^ + + p,)\diYu^f)dxds <Co. (12) 
Jo Jrs ^ ' 

• The continuity equation is satisfied in the sense of renormalized solu- 
tions, i.e.: 

dtb{p^) + div(6(/)u) + (6'(/)/ - b{p^)) divu^ = 0, 
for any b G C^(R^) such that 

b'{z) = constant, for any z large enough, say z > M. 

• The system ([U]) holds in V'{{0,T) x R-^). 

Beside the results on the existence of weak solutions for the Cauchy prob- 
lem for the Navier Stokes Poisson system (jlip there is a theory concerning 
the global existence of classical solutions of (|lip see for example |21j for the 
iJ^ framework or for global solutions in Besov spaces. We describe this 
global existence result in the following theorem. 
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Theorem 3.2. Assume that (p^ - l,mo) G H^{R^) H L^{R^), s > 4, with 
8 = IKpo "^o)||H''(R3)nLi(M3) small. Then, there is a unique global classical 
solution {p^,m'^,V'^) to the system (jlip satisfying 

p^-ie c°{R+,H'{M.^))nc\m+,H'-\M.^)), 



3.2 Main results 

Having collected all the preliminary material we are now ready to state our 
main results. The first result concerns the convergence of solutions of the 
system ([TT]) in the quasineutral regime. 

Main Theorem 1. Let [p^jU'^jV"^) be a sequence of weak solutions in 
of the system ill]), assume that the initial data satisfy (ID). Then 



(i) / — > 1 weakly in L°°([0, T]; 

(ii) There exists u£ L°°{[0,T];L^{M.^))nL'^{[0,T];H'^{R^)) such that 

weakly in L'^{[0,T];H^{M.^)). 

(iii) The gradient component Qu^ of the vector field satisfies 

Qu^ — > strongly in L'^{[Q,T];L^{R?)), for any p e [4,6). 

(iv) The divergence free component Pu^ of the vector field satisfies 

Pu^ — > Pu = u strongly in L'^{[0,T]; Lf^^{R^)). 

(v) There exist correctors , E~ in L°°((0, T), L^(M^)) and a positive 
microlocal defect measure on R? x 5^ depending measurably on t, 
associated to the electric field E'^ = \IV^ , such that for all pseudod- 
ijjerential operators A G ip^^^plR"^ , ]C{R^)) , and of symbol a{x,^) and 
for all (j) G T'(0, t) one has 

lim^y dt(j){t)X^{AE^,E^) = j dt(j){t){AE+,E+) + J dt(j){t){AE- , E'] 

+ / dt<l>{t)[ tr(aix,C)^)du''. 

(13) 
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(vi) u = Pu satisfies the following equation 



Pidtu - Au+ {u-V)u- 



div{E+ <S)E+ + E-^ E-) - div(i/-^, ) = 0, (14) 

in P'([0,r] X M?). 

Remark 3.3. In the previous theorem we constructed a defect measure 
and the correctors E^. They correspond to the physical phenomenon of the 
high frequency plasma oscillation. Notice that the correctors E^ remain 
important as A — t- and are not vanishing, in fact we don't have initial layer 
but on the contrary the effect of ill prepared initial data appears through 
E^ and remains important for all times. 

As we will see in the rest of paper the construction of the defect measure 
will be done by using the theory developed by P. Gerard in [12] and L. 
Tartar in [29]. The explicit construction of the correctors is not trivial and 
requires a smooth setting for the solutions. We will show this part in the 
next theorem. 

Main Theorem 2. Let be {p^ ,V'^) be a sequence of the Navier Stokes 
Poisson system, satisfying for s > 4 

11/''^ — l||L°°(0,T;_ff'=(M3)) < C ||Ail^'^||z,oo(o,T;J?''(R3)) < C (15) 

then, for all s' < s — 2 
i i 

X{E^ - e-'*/^E+ - e'^/^E-) strongly in C^{0,T, H^^~^{R^)). (17) 
and E^ satisfy 

dtE^ - AE^ + Q div{v (g) E^) = 0, PE^ = 0. (18) 

In the Section 8 we will show in the Proposition 18.41 the existence of 
solutions for the equation (fT5|) . The rest of the paper is devoted to prove 
the Main Theorems [T] and [2j 



u 



A _ _e-it/A^+ _ _eiVA^- _^ ^ strongly in C^{Q, T, H^^~\R^)). (16) 



4 Uniform estimates 

In this section we wish to establish all the a priori estimates, independent 
on A, for the solutions of the system ()lip which are necessary to prove the 
Main Theorem [TJ 
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4.1 Consequences of the energy estimate 

We start by collecting all the a priori bounds that are a consequence of the 
energy inequality ()12p . Before going on let us define the density fluctuation 

as 

0-^ = / - 1. (19) 

Proposition 4.1. Let us consider the solution {p'^,u'^,V^) of the Cauchy 
problem for the system pT|) . Assume that the hypotheses (ID) hold, then it 



follows 










is bounded in L°°([0, T]; L^(II 


■l^)), where k = min(7,2), 


(20) 




is bounded in L'^{[0,T] x M^) 




(21) 




is bounded in ^^([O.r] x M^) 


nL'^{[0,T]-L^{R^)), 


(22) 




is bounded in L'^{[0,T]; ( 




(23) 




is bounded in L'^{[0,T]; 1^(1 




(24) 



Proof From ([l2|) it follows that vr^ G L°°{[0,T]; L^{R^)). By taking into 
account that the function z — )■ z''' — 1 — 7(2; — 1) is convex and by following 
the same line of arguments as in [23] we get when 7 < 2 that 

sup / 1 1/ - l|'x|pA„i|<i/2 + \p' - irX|pA„i|>i/2| (t, x)dx < C (25) 

and when 7 > 2, 

sup / \p^-lf{t,x)dx<C, (26) 

so we can conclude that is uniformly bounded in A in L°°([0, T]; L2(M^)), 
where k = min(7, 2). (|21|) and (|24p are a consequence of (|12|) . The fact 
that G L'^{[0,T];L^{R^)) follows from 1^ and by Sobolev's embeddings. 
Now we prove E L'^{[0,T] x M^). 

/ \u^\'^dx= |k^|^X|pA„i|<i/2 + |^^^|^X|p=-i|>i/2| c^^; 
<2 [ p^\u^\^dx + 2\\p' -l\\Lk\\u^\\\,/,., 

<Co + Co\\u'\\lY\\Vu^\\l, (27) 

We complete then easily the prove by using (j21|) . Recalling that 7 > 3/2 
and by interpolating we get that E L^{[0,T]- L'^{R^)nL'^^/^^~'^\R^)). By 
using ([20]) we obtain that p-^u-^ is uniformly bounded in L^([0, T]; L^/^(IR^) + 
2^2fc/(fc+i)(-]^3))_ Therefore by Sobolev's embeddings we get ([MD- □ 

We want to complete this paragraph with a remark concerning the reg- 
ularity of the initial data. 
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Remark 4.2. With the same procedure as for u^, taking into account (ID) 
we get that CTq is bounded in L^M'^) hence in ff~^(R'^), since 7 > 3/2. If 
we rewrite uiq in the following way 



A 



we get that rriQ is bounded in L2(M3) + ^2fc/(fc+i)(]g3) ^^j^^ ^^^^^^g ^-1 
Finally we can conclude that 

CTg, ttiq are bounded in i/~^(R'^) uniformly in A. (28) 
4.2 Density fluctuation acoustic equation 

From the estimates of the Proposition 14. II we get only the weak convergence 
of the velocity field and unfortunately this will be not sufficient to pass into 
the limit in the nonlinear terms (such as the convective term div(p'^ti^ 
u^)) of the system In particular this weak convergence is induced by 

the rapid time oscillation of the acoustic waves or by the so called plasma 
oscillations. In order to overcome this problem we will estimate the density 
fluctuation uniformly with respect to A. So we derive the so called 
acoustic equation which governs the time evolution of a'^. First of all we 
rewrite the system ()12p in the following way 

dta^ + div(p^n^) = (29) 
dtip^u^) + Vcr^ = fiAu^ + + Ai)V divii^ - div(p^M^ ® u^) 

- (7 - l)V7r^ + cj^VF^ + Vy^, (30) 

X^AV^ = a^. (31) 

Then, by differentiating with respect to time the equation ()29p .taking the 
divergence of ([50]) and by using ([HT|) we get that o"'^ satisfies the following 
equation 

dtta^ - Acr^ + — = - div(^A'u^ + (i^ + /i) V div u^) (32) 
A^ 

A I „Av7T/A 



+ div [dw{p-^u-^ (g) n-") + (7 - l)V7r^ + a-^W 

It turns out that (j32|) is a nonhomogeneous Klein Gordon equation with 
mass 1/A. In order to get some more uniform estimates on we apply to 
P2p the Strichartz estimates ([8]). To renormalize the mass of the equation 
p2p more easier to handle we rescale the time and space variable, the density 
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fluctuation, the velocity and the electric potential in the following way 

r=|, . = f (33) 

u{y, t) = u^{Xy, At), t) = p^{Xy, At) 

a{y,T)=a\Xy,XT), V{y,T) = V\Xy, Xr). (34) 

As a consequence of this scaling the Klein Gordon equation (|32p becomes of 
mass equal to one, namely 

drrO' — Afj + a = —— div(/iAn + {u + /_f) V div u) 
A 

+ div (div{pu u) + (7 - l)V7r + aVV^ . (35) 

Now we consider cj = cti + (J2 + (J3 where cJi , cj2 , cJa solve the following Klein 
Gordon equations 

{drrO'i — AcJi + ai = — J div(;uA{t + (z> + /i) V div u) = Fi 
(36) 
(7i(x, 0) = a{x, 0) = (To 5t-(7i(x, 0) = dr^{x, 0) = dtdo, 

{dTT^2 - A(T2 + 5-2 = div(div(p'u ® u) + (7 - l)V7f) = F2 ^^^^ 
(T2(X, 0) = (?ra-2(x, 0) = 0, 

1 a^.ag - Aag + as = - div(aV V) = F3 ^^^^ 

|<T3(x,0) = (9^(T3(2;,0) = 0, 

We are able to prove the following estimate on a^. 

Theorem 4.3. Let us consider the solutions (p^jU'^jV^) of the Cauchy 
problem for the system pip with initial data satisfying (ID). Then for any 
So > 3/2, the following estimate holds 

^"^ W^'^W Lfwr^-'-* + \\9t(r^\\Liwr°-''' 

< A^«-^||cTo^||^-3/2 + A^«-^||m^||^-5/2 

+ r|| div(div(p^u^ u^) - (7 - i)v^^)||^^^^-.o-2 

+ A'o II div Au^ + V div ^^11^2^-2 + r|| div(f7^y^) ||^^^-.o-2 . (39) 

Proof. Since o"i, a2, o"3, are solutions of the equations (fH7|) . PS|) . we 

can apply the Strichartz estimate (l8|) with (y,T) G x (0,T/A). We start 
with cji. From ([21]) we deduce that Fi £ L^H~'^, so by using ([5]) with s = 2 
we get 

II^iIIl^H^-^'* + \\9r^l\\LiWy'^* ~ 1 1^0 1 1^:^-3/2 + 115,^0 11^-5/2 

+ A"^r||A~^/2 div(Au + V div u) 11^2 j^-2. 

' ^40) 
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From the estimate (|12p we have that /jlup, vr G L'^Ll., but is continuously 
embedded in H~^°, sq > 3/2, so we have that F2 G L'f°H~^'^~^. If we apply 
dH]) to we obtain for any sq > 3/2 

ll^2|li4^y-0-2.4 + ||a^a2||^4iy-.o-3,4 

< A-i/2r|| div(div(pu ® n) + Vf^)||^^^-.o-2. (41) 



By using the Poisson equation ()3ip we can rewrite -F3 as F3 = div(div(Vy (8) 
Vy)+iV|Vy|2). Taking into account (USD, as for F2, we get F3 G L^if-*"-^ 
for any sq ^ 3/2. Hence o"3 satisfies 



< X-^/^T\\ div(VF C3 Vy) + Jv|Vl/p|L „-.o-2. (42) 



Summing up (gO]), dH]), (02]), a verifies 

W^hiWr + l|5rCT"|li4^-»0-3.4 < ll^0||j:^-3/2 + 1 1 5,^0 1 1 ^-5/2 



+ A~^||A"^/2 div(Au + Vdivi2)||^2J^^- 
+ A~-^r|| div(div(p-u (g) n) + Vvf) 



-1 



+ A-ir||div(aVy)||^^^-.„-i (43) 

Finally, since 

ll~ll + £11 

and by using (f28|) . we end up with ([39]) . □ 



5 Strong Convergence of the velocity field 

In this section we will study the strong convergence of the velocity field 
u^. This will be achieved by studying separately the convergence of the 
divergence free vector field Pu^ and of the gradient vector field Qv^ . 

5.1 Strong convergence of Qu^ 

Here we prove the convergence of Qv?" to 0. The main tool in this process 
lies on the fact that Qv?^ can be computed in terms of o"^, so we can use the 
estimate ([5^ combined with the Young type inequalities Q , (fTU|) . 

Proposition 5.1. Let us consider the solution {p^,u'^,V^) of the Cauchy 
problem for the system (jlip . Assume that the hypotheses (ID) hold. Then 
as X I 0, 

Qu^ — > strongly in ^^([0, T]-LP{W^)) for any p G [4, 6) . (44) 
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Proof. In order to prove the Proposition 15 . 1 1 we split Qu"^ as follows 

where ja is the smoothing kernel defined in Lemma 12.31 Now we estimate 
separately Ji and J2. For Ji by using Q we get 

Ji < a^~^(^"^)||Vu^||^2 . (45) 



To estimate J2 we take into account the definition ()19p and so we split J2 

as 

J2 < \\Q{a^U^)*ja\\L^L% + \\Q{p^U^)*ja\\L^L^^ = J2,l + ^2,2- (46) 

For J2,i we use ([25]) and PU]) . so we have 

J2,i < a"^"^(^"^) ||cj^u^||^2^-i. (47) 

From the identity Q{a^u^) = V A^^dtcr"^ and by the inequality (jlOp we get 
J2,2 satisfies the following estimate 



J2 = Ai/2||A-i/2vA-i9ta^*j||^2^P 



< Ai/2a-^«-'-<^-^^)rV2||A-V2a^^.|| ^ (48) 

Wx 

Now, summing up ([17|) and we get 

IIQn^llL?LS < a^-K^i) + CT\'/'a~''~'<^^^ , (49) 

Finally, we choose a in terms of A, for example in a way that the two terms 
on the right-hand side of the inequality ()49p are of the same order, namely 



a = A"+4so. (50) 

Therefore, we obtain 

IIQ^^IlLfLP < CrAf^i'+^'-o) for any p G [4, 6). 

□ 
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5.2 Strong convergence of Pu'^ 

It remains to prove the strong compactness of the incompressible compo- 
nent of the velocity field. To achieve this goal we need to recall here, The 
compactness can be obtained by looking at some time regularity properties 
of Pu^ and by using the Theorem 12.41 but before we need to prove the 
following lemma. 

Lemma 5.2. Let us consider the solution {p'^, u'^, V'^) of the Cauchy problem 
for the system ([TT]) . Assume that the hypotheses (ID) hold. Then for all 
h G (0, 1), we have 

\\Pu\t + h)-Pu\t)h2^[0,T]xR^)<CTh^^'- (51) 
Proof. Let us set = u'^{t + h) — u'^{t), we have 

\\Pu^{t + h) - Pu^{t)\\'i2 =11 dtdx{Pz^) ■ {Pz^ - Pz^ * ja) 

Jo Jr3 

+ [ [ dtdxiPz^)-iPz^*ja)=h + l2. (52) 
Jo Jr3 

By using ([9]) together with ()22p we can estimate Ii in the following way 

h < \\Pz^hiJ\Pz\t) - {Pz^*jamh2 < a||u^||^2 J|Vn^||^2^. (53) 
In order to estimate I2 we split it as follows 

l2= [ [ dtdxP{p^z^) ■ {Pz^ * ja) + [ [ dtdxP{a^z^) ■ {Pz^ * j^) 
Jo JR3 Jo Jr^ 

= h,i+h,2- (54) 
/2,2 can be estimated by taking into account (|22p . (|23p and (jliP g so we have 



-^2,2 — 



/ / dtdx{AV^z^{Pz^ * ia) 
Jo Jr? 

= f [ dtdx \{VV^z^){VPz^ * ja) + VV^Vz^{Pz^ * jo 
Jo Jrs L 

< A||AVF^z^ + AVF^Vn^||^2^i||VPz^*i„|| 

< Aa~3/2||v-u^||^2 JIAVF^z^ + XVV^Vu^W^Li- (55) 

Now we estimate /2,i. Let us reformulate P{p^z'^) in integral form by using 



18 



the equation (jlip n and the Poisson equation (jlip o. hence 



hi < 



T 



t+h 



dt I dx I ds {div {p^ (^u^) + Au^){s,x) ■ {Pz^ * j a) {t,x) 



+ 



rt+h / A \ 

dx / dsP i—VV^] {s,x) ■ {Pz^ * ja){t,x) 

fT p r-t+h 

/ dt I dx I ds{div{p^u^ ® u^) + Au^) ■ {Pz^ * ja){t, x 
Jo Jr^ Jt 



+ 



t+h 



dt dx dsX'dw{VV^0VV^)){s,x)-{Pz^*ja){t,x) 



(56) 



Then, by integrating by parts, by using ([TO]) with s = 0, p = oo, q = 2, we 
deduce 



Hi?. IIP 1^^ 



.Ai ,A|2 



iL^Ll + imV^l'hrLl)- (57) 



Summing up /i, l2,i, ^2,2 and by taking into account ()12p we have 
||Pn^(i + /i) - i^^^^(t)|li2([o,T]xR3) < C(/i2/^ + h + ha~'^/'^T^'^ + 2Q-3/2A), 

by choosing A < a and a = h?/^ , we end up with (|5ip . □ 

Corollary 5.3. Lei us consider the solution (p^, n"^, 1/"^) o/ i/ie Cauchy prob- 
lem for the system ()lip . Assume that the hypotheses (ID) hold. Then as 

Pu^ — > Pu, strongly in L^(0, T; Lf^^{M.^)). (58) 

Proof. By using the Lemma 15.21 and the Theorem 12.41 and the Proposition 
[Owe get (152]). □ 



6 Convergence of the electric field 

This section is addressed to the study of the convergence of the electric field 
= VV^. By the a priori estimate ()24p we only know that XE'^ is bounded 
in L^^LI. which does not give enough information to pass into the limit in 
the quadratic term p^VV^ = div(A^^ ® XE^) - l/2V\XE^\'^, appearmg m 
the righthand side of (jlip a. Hence the problem is how to recover the weak 
continuity of quadratic forms in . Since XE^ is bounded in LfL^. we 
can define the so called microlocal defect measure introduced by P. Gerard 
in [12] and by L. Tartar in [29j (H- measures) , but in order to handle time 
oscillations we need to introduce correctors. In this section we will be able 
to prove the following theorem. 
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Theorem 6.1. Let be {p^,u^,E'^) a sequence of solutions of the Navier 
Stokes Poisson system (jlip . then 

i) there exists E+ , E' in L°°{{{},T),L'^{M?)), 

a) there exists a positive measure on x depending measurably 
on t 

such that for all pseudodijjerential operators A G ip^^^piM.^ , }C(M.^)) , and of 
symbol a{x, ^) and for all (j) € 2^(0, t) one has 

limy dt(l){t)\^{AE^,E^) = j dt(t){t){AE^ , E+) + j dt(l){t){AE~ , E~) 



+ / dt^it) I tr ( a{x,0^^ ] du^. (59) 
First we rewrite (j32p in terms of E'^, namely 



X'^dttE^ + E^ = div A-^V div (/ti^ (g)u^ + {p^^I - X^E^ (g> E^ 
A2 



+ y div (^l^^l^lj - 2Vdivn^ = F^, (60) 

then we observe that by using (j39p and the uniqueness of the weak limit we 
have 

XVV^^O weakly in l2(0,T;L2(m3)). (gi) 

By (|6ip we see that we are exactly in the framework described by P. Gerard, 
but we have to pay attention to one fact. In our case in the quadratic 
form X'^{AE'^, E^), ^ is a pseudodifferential operator homogenous only with 
respect to the x variable and in the general case we cannot extend it to a 
pseudodifferential operator homogenous in {x,t). Hence we have to work 
on XE^ in order to isolate the components that oscillates fast in time, for 
that reason we introduce what we call the correctors of the electric field. By 
using (j60p and Duhamel's formula we can write the electric field E^ as 



£Me.</A + £Me-a/\ (62) 

where and £2 are two functions in defined by the initial data of E'^ . In 
order to understand how to isolate the oscillating terms let us consider the 
equation (j60p in the case when F'^ does not depend on x and =£^ = 0. 
Then, if we take the Fourier transform with respect to time we have {E 
denotes the Fourier transform with respect to time) 
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we can see that all the L^-mass of XE"^ is concentrated in r = itl/A as 
A —7- 0. This simple facts leads us to introduce correctors in time of order 
1/A. So we define 

El = \e-'^'^E^ E^ = Xe'^/^E^ (63) 

In particular they take into account of the L -mass of \E^ around 1/A. 
By construction it easily follows that E^ and E'^ are bounded in L^^ and 
converge weakly to E'^ and E~ respectively. Moreover we have 

Lemma 6.2. Let he {p^, u^, E^) a sequence of solutions of the Navier Stokes 
Poisson system (jlip which satisfy (ID), then one has 

£i{x) + r ds^^^^^^e^''/^ ^ E+ m P((0, T) x M^). 
The same holds for E~ . 

Proof. The proof follows by using ()62p and Proposition 14.11 □ 

So, if we look at the limit of XE^ - e^^l^E^ - e'^^I^E' as A ^ 0, we 
expect to take away the L^-mass of XE^ which concentrates around 1/A. 
Now we can define 

E^ = E^- e''/^ — - e-^*/^ — , (64) 
A A 

then we can prove the following lemma. 

Lemma 6.3. Let he {p^^u^^E^) he a sequence of solutions of the Navier 
Stokes Poisson system (jlip which satisfy (ID), then it holds 

XE^-^0 weakly in L^{0,T,L^(R^)). 

Proof. The proof follows by taking into account (|6ip and that XE'^ is bounded 
in .J,. □ 

At this point we can hope that the weak convergence of XE^ is caused 
only by spatial oscillations, which allow us to introduce the microlocal defect 
measure in space. In order to do this, since the solutions are defined only in 
(0, T), we need to extend E^ and to out of this interval and to cut-off 
the frequencies greater than a certain quantity. This will be done in the 
next proposition. 

Proposition 6.4. With the same assumption as in Lemma \6.S\ we have 

dx / dt\XTE^{t,x)\ 0, (65) 



'l?l<-R 

for any R independent on X. 
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Proof. Let be XR the characteristic function on B{0,R) and let be Tr = 
J^~^XR^ the operator that cuts the frequencies greater than R, clearly Tr 
is a bounded operator from to H^, for any s > and VTr=TrV . If we 
apply Tr to ([62]) we have 

-r pA,, , [' TrF\s,x) . ft-s\ 
TrE {t,x) = J^ 2a [-T- ) 

and by the estimates of Proposition 14.11 we have that TrF'^ , TrE^ , TrE2 
are bounded in L^L"^. Since the solutions that we consider are defined in 
the time interval < t < T, in order to use the Fourier transform in time 
we need to extend them to 0, so we get that TrF^ is bounded in L^^. Let 
us introduce 

Hl= f TRF^{s,x)e-'''^ds + 2irRE^{x). (67) 

JO 



If we compute the space and time Fourier transform of ()67p . for r large 
enough we get 

Tt,xHl{T, i) = ^Ft,xTRF^ {t + (68) 
so Ft^xH^ is in in a neighborhood of |r| = oo, moreover we have that 

/ / \Ft,xHl{T,0?dT<^ f [ \TRF\t,x)\^dxdt. (69) 

Jr3 J\t\>a ^ Jr3 

As a consequence we have that all the mass of XFi^xTrE'^ is concentrated 
in 1/A. In fact 



we have that for any rj > there exists A such that 



(70) 



/ / \XFt,xTRE^\^dCdT < T]. 



In order to proof (|65p we take into account the decomposition of E^ in 
and the following properties. For any 77 > there exists A such that for any 
A < 1 one has 

/ / \XFt,xe^''/^E^\''d^dT<r]. (71) 

JR? J\t±^\>A 
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Moreover 



I I \:Ft,,{XTRE^-e''/^E+)\dTd^ 

= 11 \Tt,.{iXTRE^e-''/^ -TRE+)fd^dT. (72) 

JM| Jt\>A 

and so for E~ . On the other hand we know that 

Xc-'^/^TrE'^ = Xc-'^I^TrE^ - TrE+ - c-^'^/^TrE-, (73) 

the same holds for E^ . By taking into account ([7T]) . ([72]) . ([73]) and Parseval's 
identity we conclude in the following way. 



/ dC [ dt\XTE^\^ = [ di I dT\XFt,^F-\RFE>^\'' 

j\i\<R J J\i\<R J 

[ di f dT\Tt,.TRE^\^ ^ 0. (74) 



□ 



Now we are ready to prove the existence of a microlocal defect measure 
for the electric field i?'^. We start by proving the L^-orthogonality of E^, 
E- and E^. 

Proposition 6.5. For any A e ■tp'^^^p{R^ , IC{R^)) and for any (/) e V{0,T) 
it holds 

lim / dt(t){t){AE^,E^) = lim / dt(t){t){AE>^ ,e'^) 

A-^O J A^O J 

+ j dt(j){t){AE'^,E+) + j dt(l){t){AE~ ,E~). 

(75) 

Proof. First of all we observe that 

lim I dt<t){t){AE+e'^'^,E-e-'*l^) = Q (76) 



A 

Then we also have 



lini y dt(t){t)X{AE^,E+e-'^'^) = (77) 

In fact if we denote by A* the adjoint operator of A we have that A*E~^ is 
bounded in and as a consequence for any > 0, there exists -B > such 
that 

/ / dt\FA*E+\^dt < r] (78) 

J\^\>B J 

Combining ([THD with ([65]) we get ([77]). □ 
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In order to prove the Theorem 16.11 and to get ()59p it remains only to 
investigate 

hm / dtX^^{t){AF^,E'^). (79) 

The sequence XE^ fits in the framework of microfocal defect measures of P. 
Gerard on [T^ but as already explained we need his prove to our sequence. 

Proposition 6.6. Let be a bounded sequence of functions of which 
converges weakly to 0, such that for every compact set C one has 



lim I d^ I dt\Fw^{t,^)\ = 0, 
A^o Jj^ J 



IK 

then, there exists a positive measure v^"'" ori M x x S^, such that for any 
A € 'i/'?omp(^^"^) ^(^^'^)) ^'^'^ of principal symbol a{x,^) and for any cp € 'D(M) 
it holds 

lim / dt(l){t){Aw^,w^) = {u'^'^{dt,dx,dO,(t>{t)aix,0)■ 
x^o I 



To prove the previous theorem we follow the same line of arguments as 
in [12], we start with the following lemma. 

Lemma 6.7. With the assumptions as in the Propositions 1 6. b\ it holds 

lirn^Sj j dt<j){t){Aw^,w^) = 0, (80) 

limSR / dt(l){t){Aw^,w^) > 0. (81) 
A^o J 

Proof. Since A is Hermitian we have 

Q{Aw^,w^) = —{{A- A*)w^,w^), 
2i 

where A - A* £ Vc"omp(I^^ so §0^ follows by using Proposition 

(i2T]) . For the real part let be 5 > 0, then a + 5 e C'^{S*n, C{R^) and 
we can extract the square root B, namely b = (5^/^ + b'. Let be B' such 
that B' = OPS{b') and B such that B = ip{6'^/'^ + B'), ip G C^{M.^), 
If = 1 on the support of a, then OPS°{B) = OPS^{f5^/'^ + ipB') = ipb e 
V'°„„p(M3,/:(M3)). So we have that 

B*B = \f\^6 + A + R, Re V'-1p(m3, /C(]R3)), 

but then we have 

/ dt(j){t){Aw^,w^) > -6yw^\\l2+^iRw^,w^). 



We end up with ()8ip by sending A to and by using again Proposition 
(O). □ 
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Prom now on the proof of Proposition 16.61 follows the same line of argu- 
ments as in [12]. So we can apply the Propositions 16.61 to the sequence AE'-^ 
end we can conclude that there exists a positive measure such that 

lim [ dt(p{t){AXE^,\E^) = {v^{dt,dx,di),<t){t)a{x,i)). 

A^O J 

If we apply the remark in Exercise 1.5 of |12] , since is a gradient and 
we are in the finite dimensional case we have that there exists a positive 
measure such that 

v^^ {dt,dx,d^) = ^i^jU^{dt,dx), 
this ends the proof of the Theorem 16.11 

7 Proof of the Main Theorem [1] 

(i) It follows from ([20]) and i^. 
(ii) It follows from (l22]l . 

(iii) It is proved in Proposition 15. 1[ 

(iv) By taking into account that we can decompose n'^ as = Pv/^ + Qu'^ 
and by using Proposition 15.11 with Corollarv 15.31 we get 

Pu^ — > u strongly in ^^([0, T]-Lf^^{M?)). 

(v) It follows from Theorem 16.11 

(vi) First of all we apply the Leray projector P to the momentum equation 
of the system ()lip , then we have 

dtP{p^u^) + Pdw{p\^ (S) u^) = fiAPu^ + Pdiv{XE^ (g) XE^). (82) 

It is a straightforward computation to pass into the limit in the terms 
dtP{p^u^) and APn^, so, for any tp G ^([0,T] X R"') we obtain 

{P{dt{p^u^) - pAu^), if) {P{dtu - pAu),ip). (83) 

In order to study the convergence of the convective term we decompose 
it in this way 

(Pdiv(/u^(g)ii^),V?) = (Pdiv((/ - l)u^ (^u^),ip) 
+ (Pdiv(u^ 0u^),(p) 
= h + h (84) 
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The term Ii goes strongly to zero. In fact it is enough to take into 
account — 1 goes weakly to zero in L^°L2, while by interpolation 
n'^ is strongly convergent in L'^L2 ■ Concerning I2 we have as A — )■ 0, 

h = (div(Pn^ » Qu^), Pip) + (div(Qn^ (g, Qu^),Pip) 

— > {div{u (g) u) , Pip) = {P div(n (gu),ip). (85) 

Finally to establish the convergence of the term P div(AVy^ (g) XVV^) 
we have to take the limit of (Pdiv(AVF'^ Cg) XW^),ip) for any (p G 
P([0, T] X M^), but we have no strong convergence for WV^ as A — )• 0, 
so we apply the Theorem 16.11 and use the microlocal defect measure 
defined in (j59|) and we have as A — )■ 0, 

(Pdiv(AVy^ ® AVy^),V9) = {XVV^ ^ XVV^,VPip) (86) 

+ [ [ VP^^^dv^dxdt 

Jo Jr^xs^ l?r 

rT 



+ [ [ VPip{E+ (g, + E- (g) E-)dxdt. 
Jo Jrs 



So, by using together ([83]) . ([HI]) . ([85]) . ([86]) we have that u satisfies the 
following equation in I?'([0, T] x 



P{dtu - Au+ {u-V)u- 



E 



div(^+ E+ + E- tg E-) - div(z/^, \—^) = 0. 

8 Equations for the correctors: proof of the Main 
Theorem [2] 

The purpose of this section is to show how to construct the correctors E~^ 
and E~ in the case of smooth solutions for the system (|lip . In particular we 
will perform the quasineutral limit in the framework of the Theorem 13.21 and 
we will end up with the proof of the Main Theorem [2j We will divide the 
proof of the Main Theorem [2] in different steps. First of all we decompose 
the electric and the velocity fields in order to single out the oscillating parts, 
then we show the existence of the correctors and finally the equations that 
they satisfy. A similar analysis has been carried out in [16] in the case of a 
periodic domain for the Vlasov Poisson system. 

For the computations we have to perform later on it is more convenient 
to rewrite the expression ()62p for E'^ in terms of its Fourier transform, 

TE\t,0 = ij-F(.,Osin ds 

+ TE\0,0cos (j\+XTE^{0,Osm Q") . (87) 



26 



8.1 Step 1: decomposition of the electric field 

In this section we decompose the electric field in a way to isolate the time 
oscillations. First of all we define the following operator that cuts off the 
oscillations in time, for any cp £ C^(0, T, i?*), s > we set 

■H^(p{t,x) = — -/ (l){a,x)da, 0^ = 1-%^. (88) 
2vrA J^ 

Then we decompose in the following way 

= g^E^ E^ = n^E^, (89) 

clearly, Ej' is the oscillatory part of E^, while E2 is its averaged part. The 
following proposition hold 

Proposition 8.1. Under the assumptions of the Theorem lA'.jjl there exist 
three vector fields (which are gradients), E^ , E2, and such that E'^ = 
E^ + E^ with 

(ii) dtW^ = El, <C andW^^Q weakly in L^, 

(Hi) ll^2llL-//r' - 

Proof. Since F is uniformly bounded in LfH^~^ we have that H/^E^ is 
bounded in Lf^H^~^ and so we get (i) and (iii). Now, if we define 

J^W^ = XJ^E^iO, sin (^) - XTE^{0, ^ cos * 



+ da ds^^^^^ sin + J'n^E^a, i)da. (90) 

we easily obtain (ii). □ 

8.2 Step 2: decomposition and limit system for the velocity 

Now we decompose the velocity field u^, in the following way 

= v^ + W^, (91) 

where is the corrector introduced in the Proposition 18. II and we can look 
at v'^ as the velocity field without its oscillatory part W'^. For we can 
prove the following result. 



27 



Proposition 8.2. Let = — W^, then for any s' < s — 2 , and 

p'^ converge in C^{0,T; H'^ (l^^))? respectively to v and 1 and there exists a 
function H, such that v satisfies, 

diYV = (92) 
dtv + vVv- Av = Vn (93) 

Proof. First of all we can observe that v'*' is bounded in LfH^~^ and dtv'^ 
is bounded in LfH^~^, so we have that 

— y V strongly in C(0, T; i^f^c"^)' ^ny s' < s. (94) 

Now we rewrite the second equation of (|lip in terms of and W'^. 

pdtv^ + p{v^ + W^) ■ V{v^ + W^) + v{p^y 

= A{v^ + W^)+V d\Y{v^ + W^) + pE^. (95) 

By the Poisson equation divi?'*' = p'^ — 1, we have that 

/ — ^ 1 strongly in C(0, T; H'-^), (96) 

from this follows (j92p . Then, by using Proposition 18 . 1 1 we have that VF^ and 
dxW^ are bounded in and converges weakly to 0, so we can pass into the 
limit in (1951) and we conclude with (1931). □ 



8.3 Step 3: Existence of the correctors 

In this section we will identify and establish the existence of the correctors. 
First of all we introduce the operator T^, for any (p € Lf^L"^ we set 

TT^ct>it,0 = e^'^^mt,C), (97) 

By construction we have that satisfies the following properties. 

(Tl) are selfadjoint, act isometrically on L'^H'^, for all s and T^Th = 
ThT^ = I, 

(T2) if 4>^ 4>, strongly in L^, then, 0, weakly in L^, 

(T3) if (/), V G LfHI,, s > d/2, T^(t>T^i; weakly in L^,,. 

In the next proposition we prove the existence of corrector for the electric 
field E'^ and the velocity field u"^. 

Proposition 8.3. There exists two functions E^ and E~ in C^{H^~^^) such 
that for all s' < s, 

(cl) \\XE^ - ThE+ - T^E-\\ ^ 0, 

loc ^ 
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(c2) \\W^ - lThE+ - \T^E-\\^,^^,.,^ ^ 0. 

Proof. We can split in two components E^ and E^, in the following way 

\ 2 2i Jq 2iX J 

and we define in a similar way E^. We can easily verify that 

\\^±^±\\Lfm-^ = \\^±\\Lfm--^ - \\^tJ^'^±E^\\^oojjs-i < c (99) 

So we can conclude that that there exists two curl free vectors E~^, E~ such 
that 

T±E^^E^ strongly in C7°if^'^-;^i, for all s' < s. (100) 
Since we know that is an isometry we have 

E^-ThE+ ^0 and - TjJ^^;" ^ strongly in C^^iJ^';-^. (101) 
By Proposition 18.11 we know that E'^ = E^ + E2 , so by using (jlOip we get 
A^i^ - T^E+ - T^E- strongly in C^H'J~l (102) 



In order to prove (c2) we use for W^, defined by ([gOjl . a decomposition 
similar to (j98|) . namely 




we define in a similar way. From ([90]) . ()103p . (jl04p we have VF^ = 
+ + ly^. Arguing as before we have that there exists and W~ 
such that 



and VF^ - T^P^"" ^ strongly in (105) 

and 

strongly in LfW'^. 

So we get that 

- ThW+ - T^W- strongly in C^H^'^^ (106) 

The last step is to identify who are W^. Taking into account that J^T^E^ = 
iTT^W^ and PHU]) and (fTUS]) we end up with 

= -lE^. (107) 

□ 
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8.4 Step 4: Equation of the correctors 

Finally in this section we finish the prove of the Main Theorem [2] and we 
will able to show the equation ()18p satisfied by the correctors. In order to 
do this we take the equation (|60|) 

X^dttE^ + = div A^^V div (^p\^ (g)u^ + {p^yi - X^E^ ® E 

+ — div (l^-^pl) - 2Vdivn-^, (108) 



2 

we substitute the decompositions obtained in the previous sections and we 
send A to 0. Let (j) £ P((0,T) x M^), then we have 

{T^iX'dttE^ + ^^), 0) = iE\ X^duThck + T^4>) 

= {T^T^E^, X^dttT^cp + T^cp) 

= {XT^E\ jiX^T^dttT^cP + </.)). (109) 

We know that XT^E'^ E'^weakly in and we can compute 

( X^TldttT^(t> + T^(l) + (l)\ 9 ^ 1 



2idtF(j) + XdttFcj). (110) 



So we have 



{Tl{X^dttE^ + E^),ct>) = {XTT^E\ J^j{X^T^duT^(l) + (/>)) 

= {XFT^E^,2idtF(t) + XdttF<t)) 

{FE+,2idtFcp) = -2i{dtE+, (Ill) 



The next term to analyze is the convective term p u = p {v + W ) ® 

{v^ + W^), it will be sufficient to analise the terms of this sort p^v^v^, 
p^vfW^, p^W^W^, i,k = I, . . . ,3. Since v^, are strongly convergent in 
L^, by using the property (T2) we get 

T^{p\iV^) weakly in L^. (112) 

Taking now into account (c2) we have 

lim^r_^(div A-^V div(p^w,^ (g, W^) = 

- i lim T_^((div A-i V div{p\^TlE^ + p\^T^E+)) = h + h (113) 

A — ^0 

For Ii we have 

h = -i lim r:^(div A-iVdiv(/u^r^£;^)) 

A— >-0 

- i lim(e-^*/^J'(div A"^Vdiv(pVe"'*^^-^^fe )) ^ (114) 

A — ^0 
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Concerning I2 we have 

/2 = -i{TT^{{dW A-^V diY{p\^T^E+),J^cf)) 



-le 



i{diY{viE+),(P). (115) 



Finally we estimate the term p^W^^W^, again we use (c2) and we obtain 



lim(T^(p^W,^Vr,^),, 



-iJ^T^[p^{T^E^ + T^E+)i{T^E~ + ThE+)k], Tc/)) 
lim(^(^+ + e-^''/^E-)i * T{e''/^E+ + 6"**/^^-)^, 



lim((^+ +e-2**/^^-)i(e'*/^S+ + e-^*/^^-)fc,(/.) = 0. 

A-i>0 



(116) 



The next term that we have to estimate is Vdivn"^, we have 



lim(r^Vdiv 



lim 

A^O L 



(T^ V div v^,(p) + (T^ V div 



-i lim{TT^{ThE+ + T^E-),F4)) 

A — ^0 

-i \i^{iiij{FE^ + e-'^'^/^FE-),F(t)) 

A — ^0 



-i(Vdiv^+, 



(117) 



By using (|96|) we get that div A~-'^V div((/3'^)'^I) converges strongly to zero. 
It remains only to estimate the two terms div A^^V div(A^i?'*' E'^) and 
div(|Ai?'^pI). Since the arguments works on a similar way we estimate only 
the first one. Let us denote by A the operator div A^^V div and by a(^) its 
principal symbol. Then we have 



{e~'^l^a{i)T{T^E+ + T^E") * F{T^E+ + T^E~),T(j)) = 0. (118) 

If we sum up dm]), ^U^, dm]), (HIS]), (HH]), dm]), (HnD, (HUD we end up 
with the following equation for the corrector E^ 



dtE+ + dW{v(g) E+) -VdivE^ = 0, PE+ = 0. 



(119) 



Moreover, by projecting the equation (|lU8p in the divergence free space and 
following the same steps as berfore we get as A goes to the following relation 
for E+ 

Pdiv{v(g)E^) = 0. (120) 

As a consequence of pi9|) and (|12U|) we get that E~^ satisfies the following 
parabolic equation 



dtE+ - AE+ + Q div{v ® E+) 



0. 



(121) 
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On order to obtain the equation satisfied by E~ we can follow step by step 
what we have done for E~^. From the previous paragraph it is clear that the 
equation ()12ip holds in the sense of distribution, in the next proposition we 
can establish a more precise result on the existence of solution for (jl21|) . in 
particular we will see that the kernel of the Leray projector P is an invariant 
subspace for the flow of the equation (|12ip . 

Proposition 8.4. Let us consider the correctors equation 

dtE^-AE^ + Qdiv{v(E)E^) = 0, (122) 

where v G L°°(0, T; //^(M^)), 3/2 < s < 2 and the initial data satisfy 

E'^iO) G L'^{W^), P^±(0) = 0. (123) 

Then the Cauchy problem (fT22|) - (fT23|) has a unique solution E^ G L°°(0, T; L^iM^)), 
such that PE^{-,t) = 0, for any t £ [0,T]. 

Proof. The proof follows by rewriting (jl22p in integral form and by using a 
standard fixed point argument (for more detail see [30j . Chapter IV, Exer- 
cises 7.8 and 7.9). □ 
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